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Introduction 


For career astronauts and long-duration missions, the high-energy heavy-ion component of galactic 
cosmic rays may cause harmful radiobiological effects. To assess these effects, an accurate theoretical 
description of the transport of high-energy nuclei through spacecraft structures is being developed 
(ref. 1). Essential inputs to this transport theory are reaction and fragmentation cross sections for 
heavy-ion collisions at cosmic ray energies. 

In previous work, an optical model for composite particle scattering (refs. 2 through 7) had been 
solved to first order in the optical potential within the eikonal approximation to predict reaction and 
fragmentation cross sections but through indirect means. Here only the elastic channel in the coherent 
approximation, which neglects contributions of scattering to excited states with rescatterings to the 
ground state, is considered. The optical theorem allows for the calculation of the total cross section 
from the elastic amplitude, and the reaction cross section is then taken as the difference between the 
total and total elastic cross sections. Fragmentation cross sections are obtained by distributing the 
reaction cross section between the various fragmentation channels (ref. 6). This method gives good 
agreement with experimental measurements for reaction cross sections, although some discrepancies 
are seen for fragmentation cross sections which may be due to neglecting the contribution of bound- 
state inelastic transitions to the reaction cross section and/or misunderstanding ablation effects. One 
advantage of this approach is that the only nuclear structure inputs needed are ground state matter 
densities which are accurately determined from parameterization of experimentally measured nuclear 
charge densities. 

The coherent approximation for the elastic channel and the distorted wave Born approximation 
(DWBA) for inelastic transitions represent lowest order solutions to coupled-channel equations derived 
in the optical model. In this work we consider methods for obtaining higher order solutions to the 
eikonal form of these coupled-channel equations. A nucleus will have an infinite number of levels for 
which nearly all excited state wave functions are poorly known but will represent necessary inputs 
for a complete solution to the coupled-channel equations. In the forward-scattering region, a limited 
number of low lying states may lead to a reasonable truncation of the infinite channel space for a light 
nucleus, although at high energies all nuclear levels will appear degenerate so that the accuracy of any 
channel truncation will depend on the strength and momentum dependence of the excited state wave 
functions. Alternatively, solutions exact through second or higher orders in the optical potential may 
be reexpressed using sum rules to incorporate pair or higher number correlation functions. Here the 
lack of knowledge of excited state wave functions is shifted to uncertainties in the correlation fimctions. 

In the next section we review the eikonal formalism (refs. 2 and 3). Following this we discuss 
several approaches for obtaining solutions to the coupled-channel scattering amplitudes. Analytic forms 
are obtained for the second-order optical potential (bordered interaction matrix), which we apply to 
calculate angular distributions and total cross sections for elastic and inelastic scattering of p on 
and ^He on Comparisons are made between the first-order and second-order solutions. 

Eikonal Coupled-Channels Formalism 

In reference 2, a coupled-channel (Schroedinger) equation is obtained from a nonrelativistic multiple- 
scattering series for two composite particles through the use of an effective potential and the impulse 
and closure approximations. The coupled-channel equation relating the entrance channel to all excited 
states of the composite target and projectile particles is 

(vx + V’n/i (x) = X) ^rvM,nV' (x) (x) (l) 

where the subscripts n and p label the eigenstates of projectile and target, respectively; k is the 
projectile momentum relative to the center of mass; is the channel wave function; x is the projectile 
position vector relative to the target; Ap and are the projectile and target constituent numbers; 
and m is the constituent mass. The coupling potentials are of the form 


^np,n'/Lt'W ”■ {9n9^\Vopt\g^fg^f) 


( 2 ) 


where g represents the internal nuclear wave function, and the effective potential is given by 


Vi)pt — 


Ap^A'f 


aj 


( 3 ) 


where taj is the two-body transition amplitude. Equation (1) represents the complete solution to the 
nucleus-nucleus scattering system within the high-energy approximations noted. 

Matrix notation is introduced by defining the wave vector, 


'®(x) = 


t^oo(x) 
V’oi(x) 
V’io(x) 
^11 (x) 


( 4 ) 


and the potential matrix 


U(x) = HR 


■Vbo,oo(x) 
Vbi,oo(x) 
^10,00 (x) 
Vn,oo(x) 


Vbo,oi(x) 
Vbi,oi(x) 
^10,01 (x) 
Vii,oi(x) 


Vbo,io(x) 
Vbi,io(x) 
^10,10 (x) 
Vii,io(x) 


( 5 ) 


where we have denoted the reduced mass factor by hr have arranged the elements of ’®(x) and 
U(x) in order of increasing levels of target, followed by projectile excitation. The coupled equations 
are then written in matrix form 


(Vx + ’^(x) = U(x) ’J’(x) 


( 6 ) 


The scattering amplitudes for all transitions, n, to n', /x', are the elements of the scattering amplitude 
matrix which is written 

f(q) = -y^ j e"*‘'/ ’‘u(x)^(x)dx (7) 

where ky is the final projectile momentum and q the momentum transfer, 

q = k-k^ (8) 

In references 2 and 3, the eikonal approximation was applied to the coupled-channel formalism. 
Here we outline the derivation of f(q) in the eikonal approximation. In this forward-scattering 
approximation, the following boundary condition is imposed on the wave vector: 


z 


lim 

— oc 


^(x) = 



( 9 ) 


2 



where —z\s the direction of the beam source and ^ is a constant vector with a unit entry at the 
entrance channel and zero elsewhere. This boundary condition implies that no particles are scattered 
backward which is justified in high-energy scattering where forward scattering dominates. The wave 
vector is assumed to be of the form 





( 10 ) 


where ^(x) is a phase matrix whose solution we seek. The boundary condition (eq. (9)) implies 

$(x) = 0 (11) 

z — oo 

Substituting equation (10) into equation (6) with the approximations, 



lU(x)| « 

(12) 

and 


|VxU(x)|« |kU(x)| 

(13) 

leads to the solution 



(14) 


where we use a cylindrical coordinate system with cylinder axis along the beam direction such that 


x = b + z (15) 


where b is the impact parameter vector. Finally, we ignore the longitudinal momentum transfer such 
that 

q-x = q-b + O(02) (16) 

where 9 is the scattering angle, which is assumed to be small. The result for the scattering amplitude 
matrix in the eikonal approximation follows as 

f(q) = ^ J (Pb (17) 

where the eikonal phase matrix is defined by 



(18) 


and (fib indicates integration over polar coordinates in the plane perpendicular to the incident direction 
of the beam. 

We note that Feshbach and Hiifner (ref. 8) have considered the eikonal approximation for coupled- 
channel equations for nucleon-nucleus scattering. Here the ansatz of equation (10) is replaced by one 
equivalent to 



(19) 


Feshbach and Hiifner found that a solution such as equation (17) can be obtained only if the commutator 


V{h,z), r dz'V{h,z') =0 
J — oo 


( 20 ) 
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is fulfilled. Through consideration of the form of the matrix elements of U(x) in the momentum 
representation (see appendix A), we expect the condition given by equation (20) to be met in the limit 
when the longitudinal momentum transfer approaches zero. 

Approximate solutions to the eikonal scattering amplitudes were discussed in reference 3. The exact 
solutions of these amplitudes have been hindered in the past by the mathematical difficulties inherent 
in identifying the elements of f(q) with the corresponding elements of x(b) because of the occurrence 
of the exponential of the matrix ix(b). To appreciate this difficulty, consider two matrices G and B 
related by 

G = eB 


We may make an expansion such that 


G = i-hB-— + — + ... 


and the difficulty in obtaining the elei^nts of G in terms of the elements of B is now obvious because 
we must construct all the powers of B and then sum an infinite series for each element of G. In a 
similar manner the expansion of exp(ix) reveals the infinite-order multiple scattering series, including 
excitation and deexcitation of projectile and target in the eikonal limit. The exponential of the phase 
matrix has effectively summed this infinite series. In the next section, we consider possible methods 
for identifying the elements of the eikonal scattering amplitude matrix. 

Methods of Solution for Eikonal Scattering Amplitudes 

In this section we compare several approaches for identifying the elements of the scattering 
amplitude matrix as given in equation (17). The Cayley- Hamilton theorem states that every square 
matrix satisfies its own characteristic equation. The characteristic equation for the eikonal phase matrix 
is given by the determinant equation 

lx-Ai|=0 (21) 

where A represents the eigenvalues of x- For x rank iV, equation (21) represents an iVth order 
polynomial in A which must also hold for x. Through this polynomial equation in % higher order 
powers of x can be expressed in terms of lower order powers which could allow for the summation of 
the infinite-order expansion of exp(zx). We do not expect the Cayley-Hamilton theorem to lead to a 
straightforward numerical procedure for summation of this series for an arbitrary form for x. 

A second approach for_summation of this infinite series is through a similarity transformation. That 
is, if the unitary matrix A diagonalizes x such that (ref. 9) 


A^xA = D(A£) (22) 

where D is a diagonal matrix with the eigenvalues A^ of x along its diagonal, then 

exp(zx) = AD A^ (23) 

and we find 

f(q) = ^ j {A(b)D A^(b) - i} (Pb (24) 

The solution for the elements of f (q) is now dependent on knowledge of the eigenvalues and eigenvectors 

of X. 

The elements of x are, in general, complex. A condition for a complex matrix to be similar to a 
diagonal matrix is that its real and imaginary parts be diagonalizable simultaneously, that is, that it 
be normal. This is equivalent to the commutation relation (ref. 10), 

(x,x^]=0 (25) 
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In terms of real-space transition densities, the elements of x niay be written (ref. 2) as 


XnAi,nV'(*^) = / dx' p^^>{x') J dy p„„/ (x + x' + y) t(c, y) ( 26 ) 

In appendix A we show that for transitions between states of discrete values of angular momentum the 
elements of \ may be expressed in terms of reduced transition form factors, which are real, as 

XjpMp.JrMr(b) = Bj^MpBj^Mr (-1)^^+^^ 

X GjpMp,JtMt{^) exp \-i{Mp + Mt){<I)i, + tt/2)] (27) 

where Jp = Jn — Jj* = -- J^/, Mp = Mn—M^r^ and Mp = —M^i. Because of the complex 

phase in equation (27), the condition for diagonalization as given by equation (25) will not be fulfilled 
for general values of Mp and Mp. 

Finally, we consider use of Sylvester’s theorem for identifying the elements of f(q). Sylvester’s 
theorem is based on the recognition that any polynomial relation that holds for a scalar variable will 
also be true for any square matrix. In particular (ref. 11) ^if we con^der Lagrange’s interpolation 
formula to hold for a polynomial P of the square matrix B, where B has N distinct eigenvalues, 
choosing the interpolation coefficients to be the eigenvalues, Sylvester’s theorem states 

P(B) = E P(A,) n ,/ _ , 1 (28) 


The extension of Sylvester’s theorem to the infinite series represented by exp(ix) depends on con- 
vergence conditions required for the exponential of a matrix (ref. 11) which we assume to be true. 
Then 

exp(® = S n ( 29 ) 

^=1 J = 1 ^ J 


The extension to the case of degenerate roots is considered in reference 11. Here if s is the degeneracy 
of the root, then 


exp(ix) = 


1 

s - 1 


d® ^ e*^ adj x 

dA®-i As(A) 


I 


A = A3 


(30) 


where 


As(A) — (A — A3 -|_i)(A — As^-2) ... (A — Atv) 


(31) 


The solution that follows from application of Sylvester’s theorem offers the advantage that only the 
eigenvalues and not the eigenvectors of x _^re needed and should lead to a straightforward numerical 
procedure for identifying the elements of f(q). From equation (30) we expect that this method will 
not be practical if a large degeneracy occurs. In the next section, we consider a bordered form for x 
where analytic methods allow us to solve for the form of the scattering amplitudes with any of the 
three methods discussed. 


Bordered Interaction Matrix 

The bordered interaction matrix neglects couplings between all excited states, which for the elastic 
channel leads to a solution that is exact through second order in the optical potential (ref. 8). In this 
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section we consider solutions to equation (17) for a bordered interaction matrix of the form 


Xopt X00,01 X00,01 xoo,ii ■ ■ ■ 

X01,00 Xopt 0 0 • • • 

X10,00 0 Xopt 0 • • . 

X= Xll,00 0 0 Xopt • • ■ (32) 


where 

Xopt = XOO.OO (33) 

Here the density of the excited nuclear medium is taken to be approximately the same as the ground 
state. The validity of this approximation should increase with mass number and has been estimated 
to introduce an error of approximately 15 percent for ^He and 2.5 percent for (ref. 8). 

In appendix B we find the characteristic equation for a matrix of rank N of the bordered form, 
from which we arrive at the following equation: 

(Xopt ■” A) = X00,n/i Xn/i,00 (34) 

(n,^)7^(0,0) 

A direct application of the Cayley-Hamilton theorem allows us to obtain analytic expressions for all 
scattering amplitudes. Define 

^ X00,n/i Xn/ 1,00 (35) 

n,/x#(0,0) 

and from equation (34) we find 

= 2xoptA + — Xopt (3b) 

All higher powers of A may be found in terms of expressions linear in A. For example, 

= [3x^pt + T2] A + 2xoptT2 - 2x^pt (37) 

= [4x^pt + 4XoptT2] A - 3x^pt + + 2x^ptT2 (38) 

and 

A® = [5Xopt + lOx^ptT^ + T4] A + 4xoptT4 - 4xlpt (39) 

By the Cayley-Hamilton theorem, equation (34) must also hold when A is replaced by x- Consider 

exp(ix) = i + ix - ^X^ - + ^X^ + + • • • (40) 

Substituting from equations (32) and (36) through (39) into equation (40) gives 

[exp(ix)]oo,00 = 1 + *Xopt - ^(Xopt + - ^(Xopt + 3XoptT^) 

+ ^(Xopt OXopt”^^) + ^(Xopt + SXopt'T'^ + 14Xopt^^) (41) 
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Rearranging terms gives 


[exp(ix)]oo,00 


^1 + iXopt 2] ^opt 0j Xopt + • • 



T2 

IT 



(42) 


[exp(ix)]oo,00 = e*^°P‘ cos T 


which we recognize as 
Similarly 

[exp(i)^)]oo,n/u ” ^X00,n^ ^Xopt “ ^Xopt “ ”” ^Xopt “ ^Xopt^ 


(43) 


+ ^Xopt + + • • - j 


(44) 


and rearranging terms gives 


[exp(ix)]oo,n/Li — ^X00,n^ ^ 2[^°pt gjXopt 3! 5! ^ 


which we recognize as 

[exp(x)]oo,nM = »X00,nAt (46) 

From equation (17), we now identify the scattering amplitudes as 

/elas(q) = ^ / T - 1} (47) 

and . 

/oo,n„(q) = ^ / ^xoo,w dH (48) 

In reference 8 analogous results have been obtained by using the method of similarity transformation. 
We note that although x normal, this method is applicable because not all eigenvectors are needed 
to obtain equations (47) and (48). From the form of the characteristic equation (see appendix B), 


(Xopt-A)^-"|(Xopt-A)^-T-^ 


= 0 


We recognize the eigenvalues of x as 


-^1 = Xopt + T 
^2 ” Xopt 


> 


(49) 


(50) 


Xopt 

~ Xopt ” T 


It is then straightforward to obtain equations (43) and (46) with Sylvester’s theorem but tedious to 
generalize to arbitrary rank AT. 

We note that the coherent and DWBA approximations are recovered from equations (47) and (48), 
respectively, in the limit of small T. Also, in general, the function T will have a dependence on the 
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azimuthal angle (fig. 1), such that a numerical integration over this angle will be required for the 
second-order solutions. 

The function T as given in equation (35) is directly related to the pair correlation function. This 
can be seen by substitution of equation (A9) into equation (35), 



oo 

X Foni<l)FnQ{q')Fo^{-q)F^Q{-q!) 

/i or n^O 

X t{q) t{q') 


From reference 12, we have the following sum rule on the form factors 

—1 1 / 1 \ 

^on(q) ^no(q) = ^oo(q) ^oo(q') + j ^oo(q + q') + C'oo(q,q) 




(51) 


(52) 


where A is the mass number of the nucleus in question and Coo(q, qO Fourier transform of the 

pair correlation function. Analytic models for Coo(q,qO under investigation. In the next section 
we consider a numerical study of long-range correlations involving partial summation of the infinite 
summation that appears in equation (51) for 

Finally, we note that the first- and second-order solutions to the eikonal coupled-channel scattering 
amplitudes were found by approximating the form of X- We expect that higher order solutions, though 
more difficult, could be found by approximating the form of higher powers of x* 


Physical Inputs 

As a numerical study, we compare the first- and second-order eikonal coupled-channel solutions for 
p on and on scattering. The 2“^ at 4.65 MeV, 0"^ at 7.66 MeV, 3“ at 9.65 MeV, and 4+ 
at 14.1 MeV excited states of are considered. An advantage of the bordered interaction matrix is 
that the eikonal phase matrix elements may be obtained through knowledge of form factors measured 
in electron scattering experiments, so that no excited state wave functions are needed as inputs. This 
would not be true for couplings between the off-diagonal elements. The charge form factors for the 
ground and first three excited states have been parameterized previously (refs. 13 and 14) in the form 

Fch^rgeiq) = Cq^) (53) 

where the parameters B, (7, d, and m are listed in table I. Table I also lists the form factor for excitation 
of the 4~*" state at 14.1 MeV of which we have parameterized to the data of reference 15. 

The matter form factors are obtained from the charge form factors in the following manner (ref. 16) 


FA{q) = 


-^charge (^) 


where Fp{q) is the proton charge form factor given by 


Fp{q) = e-^p9V6 


(54) 


(55) 


with Tp = 0.87 fm, and Fcm{q) is a center-of-mass correction of the form 
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Fcm{q) = 


(56) 



with 


al 


<r^ > -r^ 

z-2 , 3U-1) 
2 + 2 A 


(57) 


where < > is the root-mean-square radius of the nucleus. For the ground state of "^He, we use the 

parameterization of reference 17, 




(58) 


where C\, C 2 , and (I 2 are listed in table I. 

The two-body amplitude is assumed to contain only a central piece of the usual form (ref. 2) 


t{q) = 


( 59 ) 


where the energy-dependent parameters cr(e),a(e) and B{e) are taken from reference 18 and given in 
table II. 

The Gaussian forms for the form factors and two-body amplitude assumed in our calculation allow us 
to obtain analytic solutions for all eikonal phase matrix elements needed as inputs for our calculations. 
All integrals needed are of the form listed in appendix C. 

Results and Discussion 

The first- and second-order scattering solutions and experimental data for elastic and inelastic 
scattering of p on at 800 MeV (ref. 19) and 1000 MeV (ref. 20) and for '^He on at 340A MeV 
(ref, 21) are shown in figures 2 through 12. For p— elastic scattering (figs. 2 and 6), the coherent 
approximation (dashed line) and bordered matrix (solid line) results are nearly identical in the region of 
the forward peak where single scattering dominates. We include Coulomb effects only in an approximate 
way assuming a point Coulomb interaction. A more exact treatment is needed to completely fill in the 
first minimum. (See, for example, refs. 21 and 22.) Here spin effects may also be important (refs. 17 
and 23). The effect of coupling the elastic channel to low-lying excited states is seen in the second 
maximum (figs. 2 and 6) where the bordered matrix agrees well, whereas the coherent approximation 
underestimates the data both at 800 and 1000 MeV. The sensitivity to the number of channels included 
in the second-order calculations can be seen in figure 7 where the dashed line includes only the 2^ state; 
the long-dash-short-dash line, the 2"*" and O"*" states; and the solid line, the 2“^, 0"^, 3”, and 4“^ states. 
At larger angles, agreement with the data is poor. Here the validity of the eikonal approximation is 
suspect, and the momentum transfers being probed are beyond the region where the phenomenological 
fits to the form factors and two-body amplitudes are made. For the second-order solutions, the effects 
of channel truncation, including the neglect of short-range correlations, in the T-function may be more 
important at larger angles. 

Calculations of the excitation of the 2“*", O'*", and 3“ excited states of by 800 and 1000 MeV 
protons are shown in figures 3 through 5 and 8 through 10, respectively. The dashed line is the DWBA 
and the solid line is the bordered matrix (second-order) solution. The formalism for the DWBA is given 
in appendix D. For all excited states, the DWBA and bordered matrix give similar results in the region 
of the first and second maxima. We note that although the bordered matrix contains all couplings 
to second order for the elastic channel, the cascades between excited, which are neglected, should be 
considered a second-order effect for inelastic transitions. These cascades would be most important in 
the region of the second maximum. In the region of the third maximum, we do see better agreement 
for the bordered matrix solutions as compared with the DWBA for all transitions considered. 

In figures 11 and 12 we show calculations, respectively, for elastic scattering and excitation of the 
O'*" state of for 340A MeV ^He on collisions. The experimental results of reference 21 do 
not report the forward peak with the data beginning at approximately 5®. We do not include any 
correlation effects for "^He in our calculations. The importance of correlations is expected to increase 
for the lightest nuclei (ref. 24). 

In table III we list total cross sections for all channels considered for p—^^C scattering at 340, 800, 
and 1000 MeV. The total of the cross sections is calculated by the optical theorem, and the reaction 
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cross section is taken as the difference between the total and total elastic cross sections. The first- 
and second-order results are nearly the same for all channels. This agreement is expected because our 
angular distributions show almost complete agreement between the two solutions in the forward angles 
where most of the cross section occurs. In table III we also sum the cross sections for the bound-excited 
(BE) states calculated, cr(BE). We note that cr(BE) represents only a small fraction, <5 percent, of the 
total reaction cross section, which is an indication that the neglect of the bound excited states in the 
abrasion model (ref. 6) is a good approximation, although the importance of the giant dipole resonance 
state should be estimated. 

Concluding Remarks 

We have presented several approaches for solving the eikonal form of the coupled-channel scattering 
amplitudes. Analytic forms for second-order solutions have been obtained, and their relationship to 
the pair correlation function established. Improved agreement with experimental results for p—^^C 
scattering is seen for the second-order solutions in comparison with the coherent approximation and 
the distorted wave Born approximation. Second-order effects are seen to make a negligible contribution 
to total channel cross sections. Future work should involve inclusion of spin effects and the development 
of analytic models for the pair correlation function, 

NASA Langley Research Center 
Hampton, VA 23665-5225 
September 1, 1988 
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Appendix A 

Eikonal Phase Matrix Elements in Terms of Momentum Space Representations 

In this appendix, we reexpress the real space eikonal phase matrix elements in terms of the 
momentum space representations of the transition densities and two-body amplitude. This allows 
us to explicitly show the complex dependence of these matrix elements. 

If projectile and target total angular momentum transfer numbers are defined as 


Jp ^ Jn ~ Jn’ \ 

Jt = J 

and total projection angular momentum transfer numbers are defined as 


(Al) 


Mp = Mn- \ 
Mx = - M^f j 


(A2) 


we introduce the projectile and target transition form factors and respectively, given 

by 

PjpMp (x + y + ^) = J FjpMp (q) (A3) 

and 

= j (A4) 

and the momentum-space representation of the two-body amplitude given by 

= (A5) 


into equation (26) to find 


= “^^2^ Xoo ^^ ^^JpMp{~<l) FjTMTi*i)^i^) (A6) 

In the eikonal approximation, the parallel momentum transfer is assumed to be negligible. Letting 

q-x = qx-b + g||Z (A7) 

such that qj^ and b lie in the x-y plane (fig. 1) and assuming the dependence of the form factors 
and two-body amplitude approaches zero faster than exp(— zg||z) gives 

XJpMp,JtMt = ~2i^ ^2!r)3 f_ Jq ^(q) 

(A8) 

which reduces to 

XJpMp,JtMt{^) = <^^9e“*‘*'*FjpMp(-q)-^^JTMr(q)^(q) (A9) 

where qj^ has been replaced by q. Next, we introduce reduced form factors Fj{q) defined by 


(AlO) 


FjmM = K*,) Fj(q) 
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where 


(All) 


YjmK'I;) = 

and where Pjj^ represents the associated Legendre polynomials. From figure 1 we can see that 
and using the identity (ref. 25), 


PjM = \) 


0 

(J - M)!! 


( J + M = odd) ' 
( J + M = even) 


(A12) 


then defining coupling coefficients as 


BjM = 


f(J-M)! (J + M-1)!! 
( J + M)\ (J-M)!! 


(A13) 


gives 


yjM [\^<l>q) = 

Substituting equations (A 10) and (A 14) into equation (A9) yields 


(A14) 


XJpMp,JrMri^) = BJtMt / ^-i{Mp+Mr)<l>, ^(q) 

(A15) 

Then using 

q • b = g6 cos {4>i, - (f>q) (A16) 

and using the identities (ref. 25) 


/ e”^i cos(^i-,^2) = 27T(-i)” e‘”^i J„(xi 

Jo 


and 


J—n — 

where Jn, the Bessel function of order n, gives the final result as follows: 


(A17) 

(A18) 


XjpMp,JrMrW = 

(A19) 

where we have defined 

roo 

GjpMpJtMt (*^) = / 9 dq ^(Mp+My) W Pjp (?) (^) ^(«) (A20) 

The complex dependence of the eikonal phase matrix elements is now explicitly shown in equation (A 19) 
to occur in the exponential phase and in any complex dependence in t{q). We note, also, that 
equation (A 19) offers considerable computational advantage over equation (26). 
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Appendix B 

Characteristic Equation for Bordered Matrix 

In this appendix, we prove that a bordered matrix of the form 


■^> 1,1 

h,2 

^» 1,3 • 

• h,N ■ 

h,i 

d 

0 • 

• 0 

h,i 

0 

d ■ 

• 0 

1 

• • 

0 

0 

• • d . 


of rank N has a characteristic equation of the form 






N 

t=2 


= 0 


(Bl) 


(B2) 


where the eigenvalues of Bjy are A. 

Proof: To prove that equation (B2) holds, we consider the cofactor expansion of the characteristic 
determinant of Bjy and use induction with respect to rank. The result is easily seen to be true by 
inspection for small N so our proof is^omplete. 

The characteristic determinant of B^y is written 


\^n\ = |Bjv - Al| = 0 


(B3) 


A determinant can be expanded in terms of the cofactor matrices of any of its rows or columns. 
Choosing the first row of Djy to make the cofactor expansion gives 


N 


e=i 


or 


N 


(B4) 

(B5) 


e=i 

where A, ^ is the cofactor corresponding to the element di^£, and is the related determinant. If 

equation (Bl) is used, we find 

IDivI = (6i,i - A)(-l)i+ilDi,i| + 6 i,2(-1)^+2|d^ 2I + . . . + 6i,iv(-l)^''^|Di,iv| (B6) 

By assumption we must have 

iDjvl = (61,1 - A)(d - A)^-i - 6i,262,i(d - A)^-2 - ... - 6i,^6jv,i(d - A)^"2 = 0 (B7) 

such that ,5^ , , , X xiv-i 


lDi,i| = (d-A)^ 

lDi,2l = -(-l)'+V(rf-^)^“" 


(B8) 


|Di.jv| = -(-l)'^%,iv(d-A)^-2j 

in order for equation (B2) to be correct. 
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(B9) 


Consider a bordered matrix of the form given by equation (Bl) of rank N + 1. We have 
On comparison of equations (B7) and (B9), we see that to complete the proof we must show 


id' 


1,1 


= {d-X) 


N 


id' 


1,2 


= (d-A)|Di,2| 


id' 


1,3 


= (‘i-^)|Dl,3| 


id' 


1,^ 


= (d-A)|Di,ivl 


h,N+i 


= -{d-Xf-\N+l 


That _ 

|D'i,i| = (d-A)^ 

is easily seen by inspection. Consider 





^ 2,1 

0 

0 


0 






^ 3,1 

d- 

A 0 


0 






^ 4,1 

0 

d-X 


0 




|Di, 

2 I = 











h,N 

0 

0 

■ ■ ■ d 

- 

A 


and 












^ 2 , 

1 d 

-A 

0 

• • 0 



0 



^3, 

1 

0 

d-X ■ 

• • 0 



0 




1 

0 

0 

• • 0 



0 


|D'i,2| = 



• 


. . 






^N. 

,1 

0 

0 

. . d- 

A 


0 



^AT+1,1 

0 

0 

■ ■ 0 


a 

!- A 


For |Di 2 I} we make a cofactor expansion about the last column to find 


|D'i,2l = (-l)'^+'^|D'i,2l(d-A) 


(BIO) 


(Bll) 


(B12) 


(B13) 


(B14) 
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therefore, 


|D'i,2l = (d = A)lDi,2| (B15) 

A similar result that holds for |Di 3I through is easily seen by inspection. Finally, consider 

Pl,jv+ll 


|Di,7V+iI - 


Making the cofactor expansion about the last row gives 

|Di,jv+il = -(d-A)^"^6iV+l,l 



d-X 

0 

• 0 

0 

f> 3 ,l 

0 

d-X ■ 

• 0 

0 

h,i 

0 

0 

• 0 

0 

f>N,l 

0 

0 

• 0 

(d-A) 

bN+1,1 

0 

0 

• • 0 

0 


(B16) 


(B17) 


Thus our proof is complete. 
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Appendix C 

Analytic Forms for Gjm Integrals 

In this appendix, we tabulate analytic forms for the functions Gjm (qb) as defined in equation (A20) 
when Gaussians are used for form factors and the two-body amplitude. All integrals tabulated may be 
generated from the integral (ref. 26), 


where /i is a positive integer by differentiation with respect to the parameter o. Defining 

■^a = f {qb) dq 

Jo 

these integrals are given as follows: 

4 = ^ 


2a 

3 _ 

2a2 


7^ - 
Jq ~ 




j^ = l. (i ^ 

° o3 [ 2a'^32a^) 


„-62/4o 


= .-6V4a 




J} = 




2o3 


4a 


/6=— + 

^ 2o4 4a 32a2 j 


jAa 


4 =^ e-*’V4o 
8o3 


4 = 


4 = 


4 = 


62 


8a4 

[ 4a) 

63 


16o4 


63 


4a3 

i 16a 


-62 /4a 




(Cl) 

(C2) 

(C3) 

(C4) 

(C5) 

(C6) 

(C7) 

(C8) 

(C9) 

(CIO) 

(Cll) 

(C12) 

(C13) 
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Appendix D 

Distorted Wave Born Approximation 

In this appendix, we obtain the form of the distorted wave Bom approximation (DWBA) that 
follows from use of equation (A18) for the eikonal phase matrix elements. From references 2 and 3, the 
DWBA for the scattering amplitude has the following form in the eikonal approximated, optical model 
formalism: 

Jc C 

fjpMp.JTMxi^) = J pl) 

where Xopt(b) is the elastic element X00,0o(b). Substitution of equation (A19) into equation (Dl) gives 


fjpMp,JrMri^) = 


Then with 

q • b = g6 cos {<t>b — <f>q) 
and the identities (A17) and (A18), we find 

““ fJLV 

^JpAfp^Jr^Jify* (q) ~~^^BjpAip Bjj^Mj* 


J u 


(D2) 


(D3) 


(D4) 


where we have dropped an inconsequential pheise. The form of the angular distribution for a discrete 
projectile-target transition then follows as 


^ {JpMp, JtMt) = 1^1 Z! 1 £ BjpMp BjtMt 

foo 2 


(D5) 


where and P/r are the initial and final center-of-mass momenta of the projectile, respectively. Values 
for the coupling coefficients, for small values of J, are listed in table IV. 
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Table III. Total Channel Cross Sections for p on ^^C 


Cross section 

Tiab = 340 MeV 

3lab = 800 MeV 

Tiab = 1000 MeV 

Coherent 

Bordered 

Coherent 

Bordered 

Coherent 

Bordered 

(t(EL), mb 

54.1 

53.5 

92.5 

91.3 

103.5 

102.1 

cr(BE)j mb .... 

3.5 

3.5 

7.1 

7.1 

7.5 

7.5 

cr(2"'‘), mb 

2.2 

2.2 

4.1 

4.1 

4.3 

4.3 

mb 

0.3 

0.3 

0.4 

0.4 

0.4 

0.4 

<t( 3“), mb 

1.0 

1.0 

2.5 

2.5 

2.6 

2.6 

mb 

0.03 

0.03 

0.1 

0.1 

0.1 

0.1 

<r(re), mb 

220.9 

220.3 

238.3 

237.2 

223.8 

223.0 

(r{tot), mb 

275.0 

273.8 

330.8 

328.5 

327.3 

325.1 


Table IV. Values for Coupled Coefficients 


J 

M 

Bjm* 

0 

0 

1 

1 

1 

l/\/2 

2 

0 

1/2 

2 

2 


3 

1 

vm 

3 

3 


4 

0 

3/8 

4 

2 


4 

4 

V35/128 


*Bj-m = Bjm 


0 


2 < 







Figure 2. Theoretical and experimental (ref. 19) elastic angular distributions for 800 MeV p-^^C 
scattering. 



Figure 3. Theoretical and experimental (ref. 19) inelastic angular distributions for excitation of 2“*” state 
in in 800 MeV scattering. 






Figure 6. Theoretical and experimental (ref. 20) elastic angular distributions for 1040 MeV 
scattering. 



Figure 7. Effect of channel truncation in second-order calculations for 1040 MeV p-^^C elastic scattering. 





Second-order calculation 

- - - DWBA 
O Experimental data 


Figure 8. Theoretical and experimental (ref. 20) inelastic angular distributions for excitation of 2'^ state 
in in 1040 MeV p-^^C scattering. 
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DWBA 
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// X\ ^ 


\\ ^^6^0 


Figure 9. Theoretical and experimental (ref. 20) inelastic angular distributions for excitation of 0"^ state 
in in 1040 MeV scattering. 






^cm 


Figure 10. Theoretical and experimental (ref. 20) inelastic angular distributions for excitation of 3“ state 
in in 1040 MeV p— scattering. 


mb/sr 



Figure 11. Theoretical and experimental (ref. 21) elastic angular distributions for 340A MeV a-^^C 
scattering. 
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